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Background: The Importance Truncated No-Core Shell model (IT-NCSM) has recently been shown to extend 
theoretical nuclear structure calculations of p-shell nuclei to larger model (A max ) spaces. The importance trun- 
cation procedure selects only relatively few of the many basis states present in a 'large' iV max basis space, thus 
making the calculation tractable and reasonably quick to perform. Initial results indicate that the procedure 
agrees well with the NCSM in which a complete basis is constructed for a given iV max . 

Purpose: An analysis of uncertainties in IT-NCSM such as those generated from the extrapolations to the 
complete iV max space have not been fully discussed. We present a method for estimating the uncertainty when 
extrapolating to the complete iV max space and demonstrate the method by comparing extrapolated IT-NCSM to 
full NCSM calculations up to iV max = 14. Furthermore, we study the result of extrapolating IT-NCSM ground- 
state energies to AT max = 00 and compare the results to similarly extrapolated NCSM calculations. A procedure 
is formulated to assign uncertainties for AT max = 00 extrapolations. 

Method: We report on 6 Li calculations performed with the IT-NCSM and compare them to full NCSM calcula- 
tions. We employ the Entem and Machleidt chiral two-body N3LO interaction (regulated at 500 MeV/c), which 
has been modified to a phase-shift equivalent potential by the similarity renormalization group (SRC) procedure. 
We investigate the dependence of the procedure on the technique employed to extrapolate to the complete A max 
space, the harmonic oscillator energy (Ml), and investigate the dependence on the momentum-decoupling scale 
(A) used in the SRC. We also investigate the use of one or several reference states from which the truncated basis 
is constructed. 

Results: We find that the uncertainties generated from various extrapolating functions used to extrapolate to 
the complete iV max space increase as iV max increases. The extrapolation uncertainties range from a few keV for 
the smallest A max spaces to about 50 keV for the largest iV max spaces. We note that the difference between 
extrapolated IT-NCSM and NCSM ground-state energies, however, can be as large as a 100-250 keV depending 
on the chosen harmonic oscillator energy (Ml). IT-NCSM performs equally well for various SRC momentum- 
decoupling scales, A = 2.02 fm -1 and A = 1.50 fm _ . 

Conclusions: In the case of 6 Li, when using the softened chiral NN N3LO interaction, we have determined the 
difference between extrapolated A max = 00 IT-NCSM and full NCSM calculations to be about 100-300 keV. As 
Ml increases, we find that the agreement with NCSM deteriorates, indicating that the procedure used to choose 
the basis states in IT-NCSM depends on Ml. We also find that using multiple reference states leads to a better 
ground-state description than using only a single reference state. 
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I. INTRODUCTION 

Various techniques have been developed to calculate the binding energies of the lightest nuclei, for which A < 4. 
These are the techniques of Faddeev [1, 2], the Faddeev-Yakubovsky extension [3-5], the correlated Hyper-Spherical 
Harmonics [6, 7], or the No-Core Shell Model implemented with Jacobi coordinates [8]. These methods all reproduce 
the experimental binding energies for the Triton and the a-particle, provided realistic nuclear interactions are used 
such as those that contain a three-body force. However, their extension to heavier nuclei is difficult and not well 
suited for such applications. When one performs nuclear-structure calculations for heavier systems, such as those in 
the p-shell, the methods of the No-Core Shell Model (NCSM) [9-11], the Green's Function Monte-Carlo technique 
(GFMC) [12-15] or the Coupled-Cluster technique [16, 17] are easier to implement and computationally more efficient. 

All of these techniques have started to place the realm of nuclear structure on a firm ab-initio footing. With the 
increases in computational resources coupled with the many-body techniques at our disposal, we can solve the nuclear 
Hamiltonian through various methods, and are able to reproduce experimental spectra reasonably well [18]. Recently, 
the community has been able to predict the spectrum of 14 F [19], which was later verified by experimental results 
[20] . We are also forming a better understanding of the role nuclear interactions in a many-body setting from our 
many-body calculations. For instance, using the chiral interaction ([21], [22]), it has been clearly demonstrated that 
the three-nucleon forces appearing at next-to-next leading order (N2LO) are essential for nuclear saturation and are 
required to reproduce the observed experimental level ordering [18]. 

The large model spaces that are typically encountered in nuclear-structure calculations are in part due to the 
short-range repulsive behavior of the nuclear interaction. Recently, two different unitary methods have been used to 
'soften' the short-range repulsive nature of the nuclear interactions. These are the wi ow -k (sec Rcfs. [23-25]) and 
Similarity Renormalization Group (SRG) methods (see Refs. [26-28]). In the case of the SRG, a series of unitary 
transformations are performed on the bare Hamiltonian, dictated by an RG flow equation. A repeated application 
of the flow equation leads to a decoupling of high- and low-momentum components. The decoupling is characterized 
by a momentum scale, A, which controls the degree of decoupling that the SRG procedure imposes upon the bare 
interaction. Typical soft interactions in use today have a range of 1.5 < A < 2.5 fm . The use of soft interactions in 
the NCSM changes the rate of convergence, i.e., smaller model spaces are required for convergence than those required 
for the bare interaction. Furthermore, the SRG and bare interaction NCSM calculations will converge to the same 
ground-state energy, provided the SRG procedure is carried out in a unitary manner (a discussion on SRG unitarity 
for the p-shcll nuclei can be found in [27]). In the case of the NCSM, the rate of convergence is also hampered 
by the use of harmonic-oscillator (HO) single-particle states, since they have the incorrect asymptotic behavior for 
nuclear bound-state wavefunctions (but offer the best means of separating the intrinsic and center-of-mass degrees of 
freedom) . 

To extend NCSM calculations to the mid p-shell region, or perhaps to the start of the sd-shell, is a challenging task, 
especially if we want fully converged results. By this, we mean calculations that are free of any NCSM parameters, 
such as a dependence on the size of the many-body basis. In the case of the NCSM, the main difficulty encountered 
comes from the relatively quick rise in the number of many-body basis states present in the model spaces. Recently, 
there have been some attempts to overcome this problem by returning to a traditional shell model approach [29], by 
projecting the Hamiltonian onto a smaller sub-model space, or by selecting only a small subset of the large number of 
many-body basis states present [30, 31]. The latter method is called the Importance Truncated No-Core Shell Model 
(IT-NCSM) and will be the topic of this investigation. 

The NCSM is characterized by two parameters, iV max and Ml. iV max describes the number of oscillator energy 
quanta available to the A nucleons above the unperturbed ground-state configuration, which in turn defines the 
number of single-particle states that must be included to form the many-body basis. We will often refer to the model 
space and A max interchangeably. Since we work in a HO single-particle basis, our calculations are also dependent on 
the chosen Ml. Both of these parameters will be discussed in detail in Sec. II. 

The Importance Truncated No-Core Shell Model (IT-NCSM) has recently been used to extend calculations of p-shell 
nuclei to larger A max spaces, in which it seems that the calculations have converged [32]. In the importance truncation 
scheme, a small set of basis states is chosen from the full N max space, using a procedure based on multi-configurational 
perturbation theory. The required N max space to reach convergence in the p-shcll, in which the complete basis is 
constructed, may have many billions of basis states. Although there are computers, and some codes, that can handle 
such a calculation, it is generally still an unfeasible calculation to perform routinely. When the importance truncation 
procedure is used, the actual number of basis states kept for the same -/V max space is usually 10-15 million states, a 
relatively easy calculation to perform. However, the uncertainties that arise from this truncated basis have not been 
explored in depth. In this paper, we will discuss IT-NCSM calculations of 6 Li, for which we can compare our IT results 
to the complete model space up to N max = 14. We begin with the bare NN N3LO interaction [21], which we transform 
to a phase-shift equivalent form that has been evolved to momentum scales of A = 2.02 fm -1 and A = 1.50 fm -1 . 
For each interaction, we vary the chosen HO energy, Ml, employing the values Ml = 12, 16, 20, 24 MeV. The IT 
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results are extrapolated in such a way as to (hopefully) recover the complete -ZV max space, in which all basis states 
are present. We investigate the various extrapolation procedures in great depth and discuss the uncertainties on the 
various extrapolations (see Sec. IV). Finally, we also investigate the sensitivity of IT-NCSM to using one (e.g., the 
ground state) or several reference states (e.g., excited states) in constructing the truncated basis. We will postpone 
the issues that have been raised concerning the size-extensivity properties of the importance truncation procedure 
(see Refs. [33, 34]) and choose to discuss estimates on the uncertainties of the current implementations of IT-NCSM. 

The NCSM will be discussed in detail in Sec. II. In Sec. Ill, we will discuss the importance-truncation selection 
procedure and in the ensuing sections, we will discuss various implications of the method such as the extrapolations 
used to obtain the ground state energy (see Sec. IV), concentrating on making some reasonable estimates of the 
uncertainties induced in the IT-NCSM calculations. In section VI, we present our detailed investigation on the 
behavior of IT-NCSM calculations, as various parameters are varied. In Sec. VII, we present our conclusions. 

II. THE NO-CORE SHELL MODEL 

The NCSM is an ab initio nuclear-structure technique employing realistic two- and three-body potentials. It is 
similar to the configuration interaction approach, in which a single-particle basis of HO states are used to construct 
the antisymmetric many-body basis. Two key differences arise though. The NCSM uses realistic interactions, which 
provides us with a connection to QCD (when EFT potentials are used), whereas the traditional shell-model uses 
effective interactions, which are typically tuned empirically to experimental data as is done in Ref. [35]. Secondly, 
we allow all A nucleons to be active in our model space, in contrast to the traditional shell model, in which only the 
valence nucleons are active and the core nucleons remain energetically frozen. 

Our starting Hamiltonian, which is translationally invariant, is, 

1 A (p — p ) 2 A A 
h a = ' h^I / NN,ij+ ^2 Virm,ijk + • • • , (1) 

i<j i<j i<j<k 

where m is the nucleon mass and A is the number of nucleons present. The first term represents the relative kinetic 
energy, the second term represents the nucleon-nucleon (NN) interaction, and the third (fourth) term represents the 
three nucleon (NNN) (and higher-body) interaction. In this work, we arc investigating the convergence properties 
of the IT procedure and not those of the Hamiltonian, thus the details of the NNN (or higher-body) interaction are 
not required. Furthermore, employing a NNN force would make the calculations significantly more challenging, and 
would prevent a comparison between IT results and the full NCSM calculations in the largest N max spaces. 

Before the advent of soft-interactions, such as those generated by the SRG or Ui ow -k; it was necessary to generate 
an effective interaction by using the Lee-Suzuki method [36, 37] for a given 7V max space and a specific value of Ml (in 
the future we will use the notation {N max ,Ml) as a short-hand). The effective interaction significantly improves the 
convergence properties of the NCSM calculations. The softened phase-shift equivalent potentials can be used directly 
in their bare form, without the need for generating a Lee-Suzuki-type effective interaction (one should keep in mind 
that both SRG and vi ow -k interactions are still effective). In other words, we solve the Schroedingcr equation given by 
Eq. (1) directly, obtaining the eigenvalues and eigenvectors. Moreover, the two-body matrix elements for any of these 
soft interactions are identical in various N max spaces when Ml is fixed. This allows us to generate one interaction, 
the -/V max = 14 interaction at a specified Ml value, and use it for all smaller 7V max calculations as well. 

Once an interaction is chosen, the NCSM calculations depend only on iV max and Ml. Since we construct our many- 
body basis from single-particle HO states, our calculations are dependent on Ml. Typically, some initial calculations 
are performed, in which a range of Ml values are employed. The resulting ground-state energies depend on the chosen 
Ml. However, in large enough -/V max spaces this dependence is often weak for a range of Ml values, which leads to a 
particular value of Ml being chosen for the rest of the calculations. In our calculations, we will use a variety of Ml 
values to study the behavior of importance-truncation. 

The A nucleons are active in a model space, denoted as -/V max , which refers to the number of oscillator quanta 
the A nucleons may share amongst themselves, above the lowest non-interacting energy configuration allowed by the 
Pauli-exclusion principle. For example, in 6 Li, an N max = 4 model space would allow one-nucleon to occupy the 
N = 5 HO shell, the other valence nucleon would stay in the N = 1 shell, and the remaining 4 nucleons would occupy 
the N = shell. Alternatively, the two valence nucleons could occupy the N = 3 shell, whereas the remaining 4 
nucleons remain in the N = shell. 

The most significant advantage of the NCSM basis is that we can exactly separate the intrinsic states from the 
spurious ccnter-of-mass (CM) states. The factorization of the center-of-mass and intrinsic states is possible if one 
uses a complete NCSM HO basis as well as a translationally invariant Hamiltonian. In order to distinguish the 
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intrinsic states from the center-of-mass states, one adds to the translationally invariant Hamiltonian (Eq. 1) a HO 

center-of-mass term, (3Hcm- Hcm = 4£f + ^mAil, 2 R 2 , in which R = \J2t^^ i s the CM coordinate. This is 
commonly referred to as the Lawson projection term [38]. The Lawson term shifts the uninteresting CM states on 
the order of (3MI up in the energy spectrum, leaving the intrinsic states in the low-lying part of the spectrum. In 
all of our calculations we set (3 = 5. To reiterate, the separation of CM from the intrinsic states is only guaranteed 
when the complete NCSM HO basis is employed. In that case, there exists an orthogonal transformation between 
single-particle and relative-CM co-ordinates [39]. Once the basis is truncated, as is done in IT-NCSM, we no longer 
have a complete factorization of the intrinsic wavefunction from the CM wavefunction. Calculations by Roth (see 
Ref. [40]) have suggested that the CM contamination remains small, on the order of 100 keV. Roth has also shown that 
once the extrapolation is performed to the complete N max space, that {Hcm) tends to zero, as it intuitively should. 



III. THE IMPORTANCE TRUNCATION PROCEDURE 



The idea of importance truncation in nuclear structure has been adopted from its use in quantum chemistry, in 
which it has been used with quite some success [41]. At the core of the procedure lies a parameter, k, that is directly 
related to the number of many-body basis states kept in a certain N max model space. We will briefly present the 
main results needed and limit technical discussions of the procedure as much as possible. The paper by Roth [31], as 
well as the Ph.D dissertation [42], present all the technical details of the importance-truncation procedure in nuclear 
structure, and we refer interested readers to it, if they desire more details than are presented here. 

The importance-truncation procedure is based formally on multi-configurational perturbation theory [43, 44]. In 
what follows, we will briefly describe the key equations that are used and how they follow from perturbation theory. 



A. The Selection of Many-Body Basis States via the Importance Measure, k, 



Suppose that one wants to perform an iV max = 12 calculation for a given nucleus, but, due to limitations set by 
current computer architectures (or resources), the full-space calculation, in which all the basis states in the N max = 12 
space are kept, is not possible. However, let us assume that an N max = 10 calculation is possible, in which one is able 
to calculate the ground-state wavefunction, which we will denote as |vI'ref,N max =io)- As a first-order approximation 
to the N max = 12 wavefunction, we can estimate the amplitudes of the N max = 12 basis states using first-order 
perturbation theory. 



,(1) \ |, T , v , (</> i ,|W / |1' rc f,Ar max = io) 

l^«=12,IT) = l*ref,JV mM =10) + }^ Z ^v) ( 2 ) 

^ e Ar max = i2 tv eref > 8 P 

In Eq. (2), we have explicitly denoted \ip^ _ 12 IT ), as the approximate wavefunction of the full space N max = 12 
wavefunction. The \<j> v ) are the N max = 12 many-body basis states. | x I > rof,w max =io) is our previously calculated 
reference state, which in our example we assume is the ground-state wavefunction of the iV max = 10 space. The two 
terms in the denominator refer to the single-particle energy level of the corresponding label. In our implementation, 
we always take e ro f. sp to be the lowest unperturbed energy configuration of the nucleus. In 6 Li, this corresponds to 
taking e rc f,sp = 2 * Ml, since two valence nucleons occupy the N = 1 shell. We neglect the zero-point motion of the 
HO, since we only require the difference in energy of the single-particle states. Furthermore, e v = (12 + 2)Ml for 
the basis states in iV max = 12. This is a particular choice that we make and is known as the M0ller-Plesset type of 
partitioning. There are other choices that one can make for the energy-denominator, however, these do not necessarily 
have superior convergence properties over the simple M0ller-Plesset partitioning [43]. 

Note that Eq. (2) requires the matrix elements of the perturbation operator, W. A convenient definition of the W 
operator is to split the initial Hamiltonian H into two pieces, namely H = Hq + W. We define Hq to be that part of 
the Hamiltonian operator that only connects many-body basis states that lie in the space N max = — 10. In other 
words, Ho does not connect basis states from our reference space to the N max = 12 space and satisfies the eigenvalue 
equation -ffol^rof) = Crefl^ref)- The full Hamiltonian, H, does, however, connect basis states from iV max = 12 to the 
reference space. Thus, we can rewrite Eq. (2), by replacing the W operator with the full Hamiltonian, H , as follows. 



,(1) \ |, T , \ , vPv\H- M^rcf ,JV max =10/ I , v , Q \ 
Vn' =12 IT/ = l*ref,iV max =10> + > \9v) (3) 

max ' * — * F , . — F — ( ™ 
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Such a form is extremely convenient since we do not need to calculate any other matrix elements than those we 
already have to calculate for the Hamiltonian operator. Equation (3) indicates that the largest correction to the 
wavefunction is essentially determined by the amplitude of the corresponding N max = 12 basis state. The amplitude 
of the basis state is, in turn, determined by the Hamiltonian matrix element between the N max = 12 basis state and 
the reference state. This leads us to define the importance measure of a basis state, K v , to be 



l«V|#l*r 



)f,N„ 



(4) 



We can now use the importance measure as a way to set a threshold limit as to which basis states are included in 
the truncated iV max = 12 space. A typical value for k is on the order of a few 10~ 5 . If we now set the threshold value 
of the importance measure to some value, say 3 x 10 -5 , we only keep those basis states (</!>„) in N max = 12 for which 
k„ > 3 x 10 -5 . Since some states have an importance measure lower than this threshold, we will discard those states, 
and thus, start to truncate the iV max = 12 space. In reality, the number of states discarded depends not only on the 
threshold value, but also on which iV max basis space is currently being evaluated. Typically, the largest N max spaces 
are most heavily truncated, whereas in the first few -/V max spaces most basis states are kept. This observation agrees 
with our intuitive notion that the components of the ground-state are dominated by basis states found in the lower 
oscillator shells. 



B. Properties of Importance- Truncation and a Posteriori Corrections 

The selection procedure of the many-body basis states that are kept in the importance-truncation calculation is 
based on a first-order perturbation theory result (for the wavefunction in iV max = 12). Using a particular value for 
k leads to a specific number of many-body basis states being kept that span the now incomplete iV max = 12 space. 
We diagonalize the Hamiltonian, H, in this incomplete -/V max = 12 space. The diagonalization results in a ground- 
state energy, E^j., associated with the truncated wavefunction, |\&jv max =i2,K;)- Note that this wavefunction, which 
results from the diagonalization of the Hamiltonian in the truncated iV max = 12 space, is not the same wavefunction 
=12 k)j w hich is what we assume is a good approximation to the actual wavefunction as shown in Eq. (3). 
Choosing a smaller value for k will result in more basis states being kept, a different truncated wavefunction, and 
thus will also result in a different ground-state energy. From the variational principle, we know that the calculated 
ground-state energy will decrease as we decrease the threshold value for k. 

It is also possible to estimate the energy contribution of the discarded states by using second-order perturbation 
theory to determine the energy correction. The first-order energy correction vanishes, since we have defined the 
perturbation operator as W = H — Hq. 



E {1) = (* ref |W|* r ef) = £ rcf - Eref = (5) 

In order to determine a non-zero quantity for the correction to the energy we need to evaluate the second-order 
correction to the energy, by summing over all the discarded basis states, as shown below. 



P (2) (^|g|^rcf,jV max =l )| 2 , . 

E o, K = - 2^ — (6) 

c ^ c rei,sp 



v— discarded 
iV max = 12 



Using this result, we can improve on the ground-state energy calculated in the truncated space, E^l by adding 

(2) 

Eq } k to it. The resulting energy, 



< +2) (7) 
has a smaller dependence on k than Eq k does, as will be illustrated later. However, Eq 1 ^ is no longer variational, 

(2) 

as it is constructed from a quantity Eq ' k , which is not a result obtained from the diagonalization of the Hamiltonian. 
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C. The extension to excited states 



So far, we have only discussed targeting basis states in the larger iV max space from a reference state that we took 
as the ground-state in a smaller N max space. We can easily extend the basis selection procedure for excited states by 
replacing the reference state with the desired excited state. In definite terms, this means that we define a re threshold 
value for each state in which we arc interested, and evaluate all basis states in the larger iV max space for each reference 
state. The re threshold value is taken as the same numerical quantity in each case. Returning to our previous example 

where we selected the basis states for N max = 12 from an A^ max = 10 reference state, we now define rel as the 
corresponding re for each of the m (ground- and excited-) states present. Note that the energy denominator remains 
the same for all the reference states used. 



re 



_ l(^|g|^k^io) 

^ref,sp 



(8) 



The number of basis states kept per fixed value of re is larger for multiple reference states compared to the number 
kept when only one reference state is used. This is expected since the structure of higher-lying states might be quite 
different to the ground-state. The overlap between different reference states with the next-larger iV max basis states, 

jv =io)' wu l ^ e VGrv different depending on the structure of the relevant reference state. 



D. Implementation of Importance-Truncation 

The Importance- Truncation procedure has been built into the No-Core Shell Model Slater Determinant Code 
(NCSD) [45], which is a multi-processor code that has its roots in the Iowa State Many-Fermion Dynamics (MFDn) 
code [46-49]. However, the current NCSD code differs substantially from MFDn. In particular, the code implements 
a hash table to look up which Slater determinants connect to each other under the action of various operators, such as 
the Hamiltonian. The relative simplicity of NCSD allows for an easy modification of the code, so that the importance- 
truncation selection procedure could be done. We will now describe how this is done in our code, specifically for the 
reference state being the ground-state (the extension to excited states is much the same). 

The reader might have formed the impression that one needs to specify only one value of re, in order to do "the 
calculation". Although it is certainly possible to calculate a good approximation to the next larger N max space 
wavefunction by using just one value for re, it is not sufficient to determine the actual energy that the complete iV max 
space would give. One needs to perform several calculations for various re threshold values, each resulting in a specific 
ground-state energy, £o,k, so that an extrapolation to re = can be performed on the Eq iK values. The extrapolation 
to -Eo,k=Oj which will be discussed in great detail in Sec. IV, yields what we assume is the true ground-state energy 
of the next-larger NCSM iV max space. Whether or not this is the case will also be addressed in Sec. IV, in which we 
compare extrapolated IT energies to the full-space energies obtained from the ANTOINE code [50]. 

In order to calculate the series of Eq iK values that we need, an efficient algorithm was developed as follows. At 
the start of every calculation we determine the smallest value of re we would like to use. Most often, we choose the 
minimum value to be re m i n = 1.0 x 10 -5 . Next, we construct all the basis states of the next-larger -/V max space and 
save those to a master file on disk. The master file is split up according to how many processors are used for the 
IT-NCSM calculation, typically ranging from 768 to 1536 processors. Each processor reads in the list of unique basis 
states (</> y ) assigned to it and determines through one Lanczos iteration, which basis states satisfy the requirement 
that re„ > re m ; n . Those basis states that do satisfy this requirement are saved to a new file along with the calculated 
value of k v . This new file holds only a small fraction of the initial iV max many-body basis states. 

Since we now have a list of all basis states that satisfy re„ > re m i n , as well as their corresponding value for re„, we are 
in a position to perform a scries of calculations, in which we now vary re. We define a series of re values, for example, 
re = {3.0, 2.0.1.0} x 10 -5 , and begin the calculation at the largest re value. All states that now satisfy re„ > 3. Ox 10 ~ 5 
are read in from the saved file and are added to the many-body basis states already present. The resulting i?o.re=3.o 
energy is saved. The process repeats, in which we now add all the basis states that satisfy re„ > 2.0 x 10~ 5 that were 
not previously added. This procedure is repeated until we have calculated all the E$ >K for all the values of re given. 
The resulting series of Eg iK values are then used to extrapolate to £7o,k=o- 

The above procedure has been specific for one iV max space. In our calculations, we employ a bootstrapping idea 
in which we apply importance-truncation to several iV max spaces in a sequential order. This is very similar to the 
IT-NCSM(seq) technique of Roth [31]. We choose to begin with a complete N max = 4 space, from which we construct 
the truncated basis in A max = 6 using the appropriate reference state. We then perform a series of calculations, as 
described above, in order to determine enough Eq, k values for the -/V max = 6 space, so that a reasonable extrapolation 
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can be made to Eq <k —q. Once we have calculated the energy for the smallest chosen re value, we use that resulting 
wavefunction, |^ r ref,JV m ax=6,« mi n): as the reference state for evaluating the N max = 8 basis states. Besides checking all 
the basis states in 7V max = 8, we also re-open the master list of all previously discarded basis states for Af max = 6, 
and check if any of those states are now kept. This point will be discussed along with other observations on the 
importance-truncation procedure in Sec. VI. When we calculate the energy contribution from the discarded basis 
states, we rc-evaluate the contribution from all the states that are still discarded. For example, once we evaluate our 
JVmax = 8 basis states, we calculate the energy contribution of the discarded states from iV max = 8 states as well as 
those that arc still discarded in N max = 6. 



(2) _ K<?M#1*rcf,W max =6)P 



<l=- E 



v— discarded 
iV max =8 



E 



v— discarded 

iVm a x=6 



t-v ^ref,sp 

K^|g|^rcf,JV m „=6)| 2 
t-v ^rcf,sp 



(9) 



This series of truncated -/V max calculations continues until the desired N max space is reached, which in all of our 
calculations is N max = 14. To summarize, note that at the end of each series of calculations in a given -/V max space, 
the wavefunction corresponding to the smallest re value is used as the reference state for evaluating the basis states in 
the next larger iV max space, since that wavefunction is the best approximation (for the specified re) to the complete 
Amax space. For each iV max space, we re-evaluate all the basis states that have been discarded in the lower -/V max 
spaces to check if any of them now satisfy the minimum re threshold. 



E. A comparison to the IT-NCSM calculations of Roth 



The implementation of IT-NCSM, which we described above, differs slightly from the implementation of the IT- 
NCSM(seq) of Roth in Ref. [31]. Two differences arise: 1) the order in which the calculations are performed and 2) 
the use of an additional truncation on the reference wavefunction l^rrfjiVmax)- We will briefly describe the differences 
between the two implementations. For a thorough discussion we refer readers to Ref. [42]. 

We have already described how we proceed with our calculations in the previous subsection. In the case of the 
Roth calculations, one value of re is specified, followed by a calculation of the ground-state energies in sequential iV max 
spaces. Once the desired -/V max space is reached, a new calculation is performed with a different value of re. Our 
method differs in the sense that we start with the largest value of re and proceed to calculate the ground-state energies 
as a function of re in a fixed- iV max space, before moving onto the next N max space. 

The additional truncation on the reference wavefunction is often referred to as the c m j n cut, in which components 
of the reference wavefunction, roughly 10 times larger than re m j n , are kept (smaller components are discarded). This 
truncated reference wavefunction is then used as the actual reference wavefunction for evaluating basis states in the 
next-larger iV max space. The implementation that we have used in this article never uses the c m i n cut. 



IV. EXTRAPOLATING IN A GIVEN iV max SPACE 



In Sec. HID, we pointed out that several values of re are used; each resulting in a ground-state energy, Eq, k . These 
are then used to extrapolate to £t),k=o- At re = in an IT-NCSM calculation, then all basis states are kept; thus, by 
extrapolating to re = 0, we hope to recover the ground-state energy of the complete iV max calculation. In this section, 
we will carefully analyze the extrapolation procedure and make a reasonable estimate of the uncertainty produced 
simply by using various extrapolation techniques. Such an analysis is new and needs to be done. We will show that 
different conclusions can be drawn from the extrapolations, depending upon how they were performed. It is not 
surprising to expect an uncertainty to be present in IT-NCSM calculations; however, we must point out that such 
results should be interpreted with care. One instance where some care should be exercised is in the extrapolation 
of a few -/Vmax calculations to the infinite space (-/V max = oo). We will demonstrate that each calculated IT-NCSM 
JVmax ground-state energy is associated with a small but finite uncertainty. These uncertainties tend to grow as iV max 
increases. If an extrapolation to N max = oo is now performed on these iV max points, one should expect an uncertainty 
to be associated with the predicted infinite result. The uncertainty on the infinite result is influenced by the finite 
Af max IT-NCSM calculations and their respective uncertainties. This naturally leads to the question, how large is the 
associated uncertainty in the iV max = oo result? 
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Figure 1. (color online) IT-NCSM calculated ground-state energies of Li in an a) iV m a X = 14 space and b) JV max = 10 space. 
The SRG-N3LO potential with a momentum-cutoff A = 2.02 fhi -1 as well as HQ = 16 MeV was used. The top two panels show 
the difference in calculated ground-state energies, when the Roth or present implementation is used. The NCSM result is shown 
at k = 0. The lower panel c) shows three different polynomial extrapolations to i?o,re=o using the -/V max = 14 ground-state 
energies as calculated by our implementation. Note that the extrapolated values are different and are spread across a range of 
about 150 keV. 



In Fig. 1, we present a series of importance-truncated-calculated ground-state energies for 6 Li in an iV max = 14 
space, using a range of k values. The top two panels show the difference between the Roth implementation and that 
used here in the iV max = 14 and iV max = 10 space. We note that the difference in the nature of these curves can 
be understood from the underlying basis states present in the truncated space. Consider the Ar max = 14 space as an 
example; and specifically the ground-state energy at k — 7.0 x 10~ 5 . In the case of our implementation, one has kept 
all the basis states from iV max = 6 — 12 at n — 1.0 x 10~ 5 , whereas in the Roth case one has only kept the basis states 
from Njobx = 6-12 that satisfy k > 7.0 x 10~ 5 . 

We have chosen 12 values of k given by the set « = {7.0, 6.0.5.0, 4.5, 3.8, 3.4, 2.75, 2.25, 1.8, 1.5, 1.2, 1.0} x 10~ 5 
which we will refer to as k— grid points. This choice is arbitrary, although we did space the smallest k values closer 
together, since we intuitively know that the smallest k values have a larger effect on the extrapolation than the 
larger k values do. Our chosen range of k, spanning from k = 7.0 — 1.0 x 10 -5 , is also to some extent arbitrary 
(one simply needs enough points to extrapolate to the ground-state energy at K = 0). Choosing too narrow a range 
could potentially affect the extrapolations in an undesired way by not capturing the general trend of the calculated 
k — specific ground-state energies. Choosing too large a range, could bias the fitted functions towards Eq >k values that 
are associated with large values of k, This brings us to our first point: The extrapolated values will depend on the 
chosen range of k and will also be influenced by the spacing of the values of k. We will address these issues in Sec. 
IV A. 

Another look at Fig. 1 c) suggests that the extrapolated values will also depend on the chosen function that is to 
be extrapolated. In Fig. 1 c), we present three possible choices; a 3 rd , 4 th and 5 th order polynomial. The predicted 
ground-state energy in the A max = 14 space has a range of about 150 keV between the 3 rd order and 5 th order 
polynomial. A priori, there is no class of functions that should be used in the extrapolations to the full-space result. 
Furthermore, we cannot make use of the Hellman-Feynman theorem, since the Hamiltonian does not explicitly depend 
on k. In other words, there is no strict requirement that the extrapolated function should have a zero-derivative at 
K = 0. This brings us to our second point: The extrapolated values also depend on the type of function used to 
perform the extrapolation. An analysis of various functions will be presented in Sec. IV B. 

The results presented in Fig. 1 are variational, as they are calculated from a certain number of basis states, 
which in turn are determined from the importance-truncation selection procedure for a given value of k, in which 
the Hamiltonian is diagonalizcd. An alternative way to fit the importance-truncated energies is to make use of the 
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second-order corrections to the energy, Eq K , as shown in Eq. (6). Formally, we know that as k -> 0, both E ( ' R (as 
shown in Fig. 1) and Eq 1 ^ should meet at the same extrapolated point. It is, thus, also possible to do a constrained 
extrapolation of these two curves, one involving only the first-order energies E^ K , the other including the second-order 
corrections, Eq 1 ^ , in such a way that both curves meet at the same point when k = 0. Such an extrapolation is 
shown in Fig. 2, using the same NCSM parameters as in Fig. I. Although we have only shown the fit for the 5 th 
order polynomial, it should be noted that using another polynomial will lead to a different extrapolated result. As 
will be shown later, when the constrained fit is used the spread in extrapolated ground-state energies is lower than 
that suggested in Fig. 1. However, a spread in the extrapolated values does remain, and we would like to characterize 
how large that spread is. 



A. Minimizing the effect on the chosen set of k values 



In Sec. IV, we pointed out that the extrapolations to -Eo,k=o depends on the chosen set of k values. In particular, 
the extrapolation depends on the range of the set, the number of K-grid points, as well as their spacing. In this section, 
we analyze the dependence on the extrapolated ground-state energy on these quantities. Such an analysis is quite 
interesting for the following reason: A different choice of K-grid points leads to a different extrapolated value of the 
ground-state energy. Typically the range of k is similar, spanning from a minimal value of a few 10 -5 to a maximal 
value of about 20 x 10 -5 . In the larger -/V max calculations, especially for the p-shcll, it is computationally expensive 
to have k < 1.0 x 10 -5 , since the number of states grows exponentially when the value of k is decreased. 

Usually one fits the 12 points, shown in Fig. I and Fig. 2, by using some specified low-order polynomial [31]. 
This, however, leads to one value of £/o,k=Oj without the ability to determine any uncertainty that is due solely to the 
extrapolation itself. One can have a feci for the uncertainty in the extrapolation by dropping the ground-state energy 
associated with the smallest value of k and then re-fitting the remaining points. The difference between these two 
extrapolations usually gives an initial estimate of the uncertainty 

We would like to provide an improved method for determining the uncertainty of IT-NCSM extrapolations. As a 
first estimate of the uncertainty produced by varying ranges and spacings of the grid-points, as well as the number 
of grid points, wc use the following procedure. We begin by choosing all possible combinations of 7 out of our 12 
available points, and for each of these (l?) = 792 sets we fit an extrapolating function to the data set and determine 
the extrapolated ground-state energy, i?o, re =o- An example of the distribution of extrapolated energies using a cubic 
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Figure 3. (color online) The normalized distribution of extrapolated ground-state energies of Li in the iV max = 14 space using 
hfl — 16 MeV. We use the chiral NN N3LO interaction, softened by SRG to A = 2.02 fm _1 . The extrapolations are done using 
a cubic polynomial fitted only to the first-order energies, as was done in Fig. 1. The extrapolated values are binned by 20 keV. 
We determine the median to be Eo,n,=o = —31.902 MeV and the standard deviation to be °"(i 2 ) — 36 keV as shown by the 
horizontal blue line. The green vertical line indicates the NCSM result. 



polynomial fitted to the first-order energies, Eq ' k , is shown in Fig. 3. After calculating all the extrapolated ground- 
state energies that result from the 792 combinations of grid-points, we bin the results in 20 keV bins. From the 
distribution we calculate the median as well as the standard deviation (indicated by the blue horizontal line in Fig. 
3). We chose the median, instead of the average, as it is a statistical quantity that is not sensitive to outliers in the 

(I 2 ) 

distribution. In the case of Fig. 3, we determine the median to be Eq ' k=0 = —31.902 MeV. The standard deviation 
of the distribution shown in Fig. 3 is c^ 2 ) = 36 keV. Note that the area enclosed by one-standard deviation below 

and above the median value in the distribution, shown in Fig. 3, encompasses roughly two-thirds of the extrapolated 
energies. 

The preceding paragraph lays the foundations of our uncertainty analysis. We repeat the above procedure for 

8 ), ( 9 J and finally ( 10 J combinations of K-grid points. For each data 
set we determine the median of the distribution, as shown in Fig. 3, as well as the standard deviation. The median for 
each data set varies by at most a few keV, whereas the standard deviation decreases as the number of combinations 
of K-grid points decreases. We determine our extrapolated ground-state energy, as well as the associated uncertainty 
in the extrapolation from the grid points by averaging over the calculated medians and standard deviations of the 4 
data sets, using 



Eq, k = 



E 



(V 2 ) 



0,K = 



E, 



0,K=0 



E, 



a 2 ) 



E 



11) 



0,K, = 



■ O 111 



) +ff ( 1 2 ) +<T O 



(10) 

(11) 



We should point out that our determination of the uncertainty which is generated from various combinations of 
k — grid points with which we associate a, is only a first attempt at determining the potential uncertainty of the 
extrapolations. The calculated standard deviation will, in general, differ if fewer (or more) data sets are used in Eq. 
(11). The important point that we want to make is that, although the uncertainty might change depending on the 
number of data sets used, the order of magnitude of the uncertainty, whether it be a few or tens of keV's will not 
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iVmax 


NCSM 


IT-NCSM 


Cubic [MeV] 


a [keV] Quartic [MeV] 


a [keV] 5 th 


-order [MeV] 


a [keV] 


NCSM [MeV] 


6 


0.198 


0.162 


-28.601 


» 


-28.601 


« 


-28.601 


w 


-28.602 


8 


1.579 


1.077 


-30.216 


2 


-30.211 


2 


-30.208 


1 


-30.213 


10 


9.693 


3.291 


-31.207 


2 


-31.204 


3 


-31.197 


4 


-31.176 


12 


48.888 


6.487 


-31.714 


10 


-31.744 


6 


-31.741 


21 


-31.713 


14 


211.286 


9.544 


-31.899 


25 


-31.971 


33 


-32.046 


29 


-31.977 



Table I. The extrapolated ground-state energies in various -/V max spaces for 6 Li (A = 2.02 fm 1 ) for Ml = 16 MeV, when using 
only the first-order IT-NCSM calculated points, Eq k . The table displays the various mean extrapolated values of the ground- 
state energy as well as the calculated standard deviations of the fits, which are indicated to the right of the corresponding 
extrapolated energy. The NCSM result, in which all basis states are kept, is shown in the right-most column. The dimension 
(in millions) of basis states are shown in the complete iV max space as well as in the importance-truncated space (column 2 and 
3). Note that the basis in IT-NCSM is drastically reduced in the larger iV max values. 



■^V max 


Cubic [MeV] 


a [keV] Quartic [MeV] 


a [keV] 5 th 


order [MeV] a 


[keV] 


NCSM [MeV] 


6 


-28.601 


» 


-28.602 


1 


-28.602 


2 


-28.602 


8 


-30.217 


2 


-30.211 


2 


-30.208 


1 


-30.213 


10 


-31.194 


1 


-31.196 


1 


-31.195 


2 


-31.176 


12 


-31.685 


6 


-31.702 


5 


-31.712 


2 


-31.713 


14 


-31.902 


9 


-31.925 


12 


-31.952 


13 


-31.977 



Table II. The extrapolated ground-state energies in various Vmax spaces for 6 Li (A = 2.02 fm 1 ) for Ml = 16 MeV, using 
the second-order IT-NCSM calculated points, E^ 2 ^ . The table displays the various mean extrapolated values of the ground- 
state energy as well as the calculated standard deviations of the fits, which are indicated to the right of the corresponding 
extrapolated energy. The NCSM result, in which all basis states are kept, is shown in the right-most column. 



change. Such an estimate does have implications, when extrapolations to N max = oo are performed. We also note 
that the standard deviation is generally smaller, when one uses the constrained extrapolations, as is shown in Fig. 2. 



B. Polynomial extrapolating functions 

In Section IV A, we focused our attention on determining the uncertainty that is generated from various combina- 
tions of n— grid points. The objective of that section was to average over many different possible choices of grid-point 
configurations. In this section, we address the choice of extrapolation function. As we had mentioned in the introduc- 
tion to Section IV, there is no a priori justification to using one function over another. One simply goes by whether 
the chosen function, once fitted to the calculated E 0Kj: lies on top of the data or not. As can be seen from Fig. 1, 
various ground-state energies are predicted, depending upon which function was chosen for the extrapolation. We 
will, thus, investigate various options that one might consider in fitting IT-NCSM calculated energies. We will use 
three different polynomials, a cubic, quartic, as well as a 5 t,l -order polynomial. 

For each selected function, we repeat the procedure outlined in Sec. IV A. Besides fitting the first-order results, 
Eq 1 }., as in Fig. 1, we also repeat the extrapolations using the same function for both the first- and second-order 

results, Eq 1 }. and Eq~£ 2 \ respectively, as shown in Fig. 2. Note that the constrained fits lead to a smaller standard 
deviation in the extrapolated ground-state energy. 

C. Estimates of extrapolation uncertainties 

We will now present our calculated uncertainty estimates on the extrapolated ground-state energy in 6 Li for the 
model spaces iV max = 6 — 14. Recall that the oscillator value is Ml = 16 MeV. The extrapolated ground-state 
energies are calculated as well as the standard deviation, which we associate with the uncertainty generated from the 
extrapolation, as explained in Section IV A. We also present the results for using various extrapolating functions. The 
extrapolated results are compared to the NCSM ground-state energies, as shown in Table I and II. 

We observe the following trends, presented by Table I, Table II and Fig. 4. The extrapolations to the NCSM 
ground-state energy for 7V max = 6 — 8 are very good, being within one keV of the NCSM result and independent 
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Figure 4. (color online) The plot shows the extrapolated ground-state energy relative to the NCSM ground-state energy (solid 
line), as well as the uncertainty, a, that we determine from variations in the K-grid points. The points to the left are the 
extrapolations that are generated from fitting just the first-order set of data, -Eq 1 ^. Those to the right are the extrapolations 
when fitting the second-order corrections, E^ 2 \ The oscillator value is hQ. — 16 MeV. 



of the function or method used. The agreement with the NCSM result is not surprising as most of the many-body 
basis states are kept in those N max spaces. Next, we observe that the uncertainty, a, increases as iV max increases, but 
that it is smaller in the larger iV max spaces for the constrained second-order fits {E^ 2 ^ ) than the uncertainty for the 
corresponding first-order fits. Note that the uncertainty stated here is from variations in the combinations of K-grid 
points. This result is also expected as now many basis states are discarded for the iV max = 12—14 spaces. However, 
note that for a given N max space that the uncertainties associated with each extrapolating function are roughly the 
same. This indicates that at least at some level, choosing one function over another, does not necessarily decrease the 
uncertainty from variations in the K-grid points. 

We also note that the mean extrapolated ground-state energy for a given iV max space can be quite different for various 
functions, varying as much as 50 keV for the N max — 14 space, when either the first- (Eq k ) or constrained second-order 

(E^ 1 ^ 2 ') results are fitted. The spread of extrapolated ground-state energies among the various chosen functions is 
usually quite a bit larger than the uncertainty associated with the variations in the K-grid points. Realistically, one 
does not have the NCSM calculations on hand, otherwise there would be no need for IT-NCSM, thus, characterizing 
the spread of the mean extrapolated ground-state energy is conceptually challenging. To illustrate this point, consider 
the results for N max = 12 in Table I. The cubic polynomial extrapolates to the NCSM result to within a keV, yet the 
quartic and 5*' i -order polynomial overestimate the result by about 30 keV, which is twice as large as the uncertainty 
from the variations in the K-grid points. If one does not know the NCSM result, one cannot make a reasonable guess 
as to which functional extrapolation is the correct one to use. Furthermore, it should be clear from the tables that 
the uncertainty from variations in the K-grid points are smaller than the spread associated with the use of different 
extrapolating functions. 



D. Asymptotically-correct extrapolating functions 

The observant reader will notice that the polynomial fits do not necessarily have the correct asymptotic behavior at 
k — y oo. The fitting routine might determine the leading coefficient of the polynomial function to be negative which in 
turn would cause the extrapolating function to decrease at a large value of k. One might propose several extrapolating 
functions that have the correct asymptotic for large k as a way to boost confidence in the fitting procedure. Examples 
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of such fitting functions are 



/(«) 
/(«) 
/(«) 
/(«) 



a + cexp(-6K 2 ) (l + g?k + fn 2 ) 
a + c (1 - tanh(6«)) (l + dn + fn 2 
a + c (1 - erf (6k)) (1 + dn + f k 2 ) 
a + ccxp(— 6k) (l + dn + /k 2 ) , 



) 



(12) 



in which b,c,d and / are treated as fit-parameters. We fix a to the ground-state energy determined in the previous 
Nmax space. For example, when we fit the points for the iV max = 14 space, wc will fix a to the ground-state energy 
calculated at k = 1.0 X 10~ 5 in the N max = 12 space. Such a constraint is possible in our calculations. There are both 
advantages and disadvantages to using these functions. The advantage is that we have fewer fit parameters than the 
quartic or 5 th order polynomials and furthermore, we have a definite constraint on the asymptotic value of the energy 
for large k. The disadvantage is that the extrapolating functions, shown in Eq. (12), have non- linear fit parameters. 
It is also true that non-linear fit parameters arc difficult to optimize globally and have the potential of leading to 
unstable fits, in which small changes in the starting values of the parameters lead to different final parameters being 
determined by the fitting routines. In our extrapolations, we typically use the parameters from previously determined 
(successful) extrapolations, in order to minimize the chance of an unstable fit. 

We repeat the extrapolations that were performed in Fig. 4, using only the first-order calculated ground-state 
energies, E^\. The result of the extrapolations is shown in Fig. 5 for the extrapolating functions shown in Eq. 12 
and should be compared to Fig. 4. We note that in some cases the calculated uncertainty (a) is quite large, notably 
at -/V max = 8. At iV max = 8 the calculated energies, Eq k , form a convex function in n. The extrapolating functions 
we have listed in Eq. 12 are appropriate for concave functions as is typically the case for the calculated ground-state 
energies at jV max > 10 (see Fig. 1). For the larger N max spaces these functions seem to have a smaller uncertainty 
than the polynomial functions used earlier. This is particularly true for the -/V max = 14 case in which the overall 
extent of the uncertainty is roughly 150 keV; the polynomial functions have an overall extent of about 200 keV. 

In the remaining part of the paper we will only use the polynomial functions to extrapolate to the ground-state 
energy in a given N max space. This choice was made out of simplicity of fitting polynomial functions. Note that no 
inherent difficulties are present in fitting the asymptotically-correct extrapolation functions as shown in Eq. 12. 



Often the final procedure in any NCSM calculation is to extrapolate the ground-state energies as a function of iV max 
to -/Vmax = oo. The purpose of this procedure is to remove the model parameters (iV max , HQ) from the calculations. The 
extrapolation to N max = oo removes the iV max dependence. The dependence on Ml is removed by the extrapolation to 
-/Vmax = oo, because, in principle, when the complete basis is recovered, no Ml dependence should remain. One usually 
chooses the HO frequency near the variational minimum of the ground-state (as we have done). This procedure was 
first used in [51], in which the bare N3LO interaction was used to determine the ground-state energy of 6 Li. We do, 
however, make the reader aware that a more satisfying extrapolation procedure is offered by considering concepts 
from EFT-theories (see [52, 53] for an EFT-inspired approach to extrapolating NCSM ground-state energies). 

In this work we will use the routinely used exponential extrapolation, in which the extrapolation to -ZV max = oo is 
done by fitting the ground-state energy as a function of -/V max to an exponential decay of the form a exp(— WV max ) + c. 
The constant c represents the ground-state energy at N max = oo. Furthermore, note that b represents a fit-parameter 
and is not related to the oscillator-length in any way. In order to extrapolate the IT-NCSM ground-state energies to 
-/Vmax = oo one needs to consider a few points. These are the following: 1) the predicted ground-state energy various 
depending on the extrapolating function used, 2) there is an uncertainty in each of the extrapolating functions, 3) 
and there is no guarantee in general that any one function is better than any other function (for example, one can 
not choose the cubic polynomial extrapolations over the fifth-order extrapolations). To elaborate on the third point, 
it is also not clear that the fit involving say 10 k— grid points should be better than the fit that used only 8 points 
(for a given function). 

We thus propose the following strategy in determining the uncertainty for the -/V max = oo extrapolations. At each 
JY max value we will randomly select one of the three polynomial functions that we have used as well as randomly 
select the number of k— grid points used in determining the extrapolated ground-state energy for that specific N max 
space. For example, we might randomly select the quartic-extrapolation that used 9 n— grid points at A max = 8, 
then select the cubic-extrapolation that used 7 ft— grid points for JV mnx = 10, and so on. This random selection is 
done for each N max value between N max = 8 — 14. Once we have selected these four ground-state energies, now as 
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Figure 5. (color online) The plot shows the extrapolated ground-state energy relative to the NCSM ground-state energy (solid 
line), as well as the uncertainty, a, that we determine from variations in the K-grid points. We have only fitted the first-order 
energies, Eq k , using the non- linear functions given in Eq. 12. The results show that the difference from the NCSM result as 
well as the spread in the extrapolated ground-state energy (a) is about the same as for the polynomial functions, shown in Fig. 
4. 



a function of iV max , we extrapolate to -/V max = oo using the exponential decay function as described in the previous 
paragraph. The procedure is done for an ensemble of 10 000 randomly selected points; increasing this number does 
not change the ground-state energy or the uncertainty at iV max = oo. Fortunately, the resulting 10 000 extrapolations 
to iV max = oo form a peaked distribution as shown in Fig. 6. In order to determine the final value of the ground-state 
energy at 7V max = oo, we simply determine the median of the distribution as well as the standard deviation, as we 
have done before. Recall that we choose the median instead of the average as the median is insensitive to outliers in 
the distribution. 



VI. OBSERVATIONS ON THE IMPORTANCE TRUNCATION PROCEDURE 

In the previous section, we addressed the fundamentals of IT-NCSM calculations. In particular, we addressed the 
choice of the extrapolating function as well as the variation in K-grid points. Here we address further questions, such 
as the dependence on the HO energy (HQ), the SRG momentum-cutoff scale (A), as well as the dependence of the 
IT-NCSM ground-state energies on the number of reference states used (i.e., targeting excited states). These are 
discussed in Section VI A, VI B and VI D, respectively. 

A. The dependence on HQ 

All NCSM calculations have a dependence on the chosen HO energy, M2, even when bare interactions are used. 
However, the dependence on Ml can be minimized for a range of values. In practice, one typically chooses an Ml 
range resulting in the lowest ground-state energy of the largest iV max space employed. In Fig. 7, we plot the NCSM 
ground-state energies as a function of 7V max for various HO energies (all basis states are kept). As N max increases, 
the dependence on Ml decreases, leading to a range of possible values of the HO energy. 

Having determined the NCSM results, we can now determine if the IT-NCSM extrapolated results show any 
dependence on Ml. In other words, regardless of the Ml value used, do we extrapolate to the corresponding NCSM 
result, or is there a systematic difference as a function of Mil In Fig. 8, we plot the difference of the extrapolated 
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Figure 6. (color online) The figure shows the distribution (in red) of extrapolated ground-state energies (Eq 1 }._ ) at -/V max = oo 
when an ensemble of 10 000 points is used as described in the text. The green distribution to the right represents the distribution 
that arises when only the cubic-polynomial extrapolated ground-state energies are considered. The blue arrow represents the 
uncertainty (a) from the median value of the overall distribution whereas the black arrow represents the full NCSM extrapolated 
value. (Ml = 16 MeV.) 
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Figure 7. (color online) The HO energy (Ml) dependence of NCSM ground-state energies for 6 Li. The momentum-decoupling 
scale is A = 2.02 fm _1 . Note that Ml = 16 — 20 MeV corresponds to the optimal HO energy for this interaction. The solid 
lines are meant to guide the eye. 
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Ml [MeV] 


12 a [keV] 


16 a [keV] 


20 a [keV] 


24 a [keV] 


Eg [MeV] 
E^+V [MeV] 


-32.512 205 
-32.364 59 


-32.244 148 
-32.162 50 


-32.041 124 
-31.963 48 


-31.881 132 
-31.777 60 


NCSM [MeV] 


-32.568 


-32.304 


-32.202 


-32.140 



Table III. The table shows the extrapolated ground-state energy at JV max = oo as a function of the harmonic oscillator energy 
HQ.. The extrapolations as well as the uncertainties are determined by the procedure outlined in Sec. V. We show both the 
first- as well as second-order extrapolations; these are the rows labeled E^ K and Eq 1 ^ 2 ' , respectively. The comparison to a 
NCSM extrapolation can be graphically seen in the lower panel of Fig. 8. 
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Figure 8. (color online) The figure shows the relative difference of the IT-NCSM extrapolated energies, E^_ (left) and Eg 1 ^^ 
(right), to the NCSM ground-state energy (horizontal curve) as a function of the HO energy, HQ. We also show the dependence 
on the various extrapolating functions (7V ma x = 12 — 14), indicating the uncertainty in the extrapolation technique. Note 
that there is a systematic drift away from the NCSM result, as Ml increases. The discrepancy increases as JV max increases. 
The solid black lines for the lowest panel indicate the uncertainty in the extrapolation to iV max = oo as determined by the 
ensemble-averaged procedure of Sec. V. 



IT-NCSM ground-state energies, relative to the NCSM ground-state energy. Various extrapolating functions (cubic-, 
quartic- or 5 th -order polynomials) are employed, using either the first- or second-order IT-NCSM extrapolated results, 
E^\ =0 or Eq 1 ^^, respectively. From Fig. 8, we can see that there is a systematic drift away from the NCSM result 
as Ml increases. The discrepancy also increases as -/V ma x increases, averaging about 200 keV from the exact result for 
the N max — oo extrapolations. In Section VI C we give some possible explanations for this type of behavior. 



B. The dependence on the SRG momentum-decoupling scale (A) 



We now investigate the dependence on the SRG momentum-decoupling scale, A. The NN chiral EFT N3LO 
potential, evolved to a momentum-decoupling scale of A = 1.5 fm -1 , has recently been used in the NCSM/RGM 
calculation of d — a scattering [54]. At this value of A, the off-shell characteristics of the two-body potential have been 
changed in such a way as to have the effect of producing a binding energy for 6 Li similar to that obtained including 
the chiral EFT NNN N2LO potential (the two-body terms are retained up to N3LO). This behavior is due to the 
non-unitarity of the SRG procedure, when only two-body terms are kept in the RG evolution. The effect for 6 Li is 
demonstrated in [27] (see Fig. 11). It is, therefore, of interest to compare two different SRG evolved potentials, in 
order to see if the importance truncation selection procedure behaves differently in the two cases. 
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Figure 9. (color online) The figure shows the relative difference of the IT-NCSM extrapolated energies, E^l =0 (left) and 

Eq^Zq (right), to the NCSM ground-state energy (horizontal curve), as a function of the HO energy, HQ,. We also show the 
dependence on the various extrapolating functions, indicating the uncertainty in the extrapolation technique. Note that there 
is a systematic drift away from the NCSM result as increases. A similar observation was made in Fig. 8. Note that we have 
used the same scale as in Fig. 8, so that the two figures may easily be compared. 



In Fig. 9, we plot the extrapolated IT-NCSM results as a function of the HO energy, Ml. This figure should be 
compared to Fig. 8 where we set A = 2.02 fm -1 . Note that the trends are similar, but the relative uncertainty to the 
NCSM results is a bit smaller for A = 1.5 fm - . The lower- A interactions are much softer, therefore, the convergence 
in iV max is much quicker. In this case, the IT-NCSM procedure selects fewer basis states for A = 1.5 fm -1 than it does 
for A = 2.02 fm" 1 , since for the softer potentials, fewer high-lying 7V max basis states are required to reach convergence. 
This point is clearly illustrated in Fig. 10 where we plot the number of basis states kept, as a function of iV max for 
both types of SRG evolved potentials. 



C. Further comments on the Kl dependence 



Figure 10 shows another interesting trend; fewer basis states are kept as Ml increases. In particular, note that the 
basis dimension is only about | of the size for Ml = 24 MeV than it is for Ml = 12 MeV. 

In Fig. 8 and Fig. 9, we showed that the IT-NCSM extrapolated results shift away from the NCSM results as Ml 
increases. The explanation lies in the definition of the importance measure k. Recall that n is inversely proportional 
to Ml. Thus, for M7=24 MeV, the matrix elements \(cf> v \H l^ref)! would have to be twice as large as they are for Ml=12 
MeV, in order for \<j) u ) to be kept as a basis state. The matrix element itself also depends on Ml : but this dependence 
must be weaker than the linear dependence in the denominator of k (since fewer states are kept as N max increases). 
Using this argument, it is not too surpizing that the IT-NCSM results would have a dependence on the HO energy 
and that, in general, the IT-NCSM results would be less reliable for larger values of Ml (if the same minimum value 
of k is used), than for smaller values. We should also point out that since the IT-NCSM basis is an incomplete iV max 
space (i.e., truncated), we no longer have a complete decomposition of center-of-mass from intrinsic states. Thus, a 
small amount of center of mass contamination is expected, which would increase as Ml increases. This issue has been 
addressed in [40]. 
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Figure 10. (color online) The number of basis states kept, as a function of Nmax, for the A = 2.02 fm -1 (top) and A = 1.5 fm" 1 
(bottom) SRG evolved N3LO interaction. Note that for the softer interaction, A = 1.5 fm~ , IT-NCSM keeps fewer basis states 
from the larger iV max spaces. We also note that as HQ increases, fewer basis states are kept for the larger N max spaces. 



D. Using multiple reference states 

One final feature we would like to investigate is the behavior of the excited-state spectrum in IT-NCSM calculations. 
In order to reliably calculate the excited states of an IT-NCSM calculation, one needs to employ a reference state for 
each state that is to be calculated. In our case, we desire to calculate the ground-state (J* = 1 + ) and the first two 
excited states, corresponding to = 3 + , + . We use as initial reference states from iV max = 4 each of the three states 
to generate the basis states that are kept in N max = 6. Since we are now using several reference states, the basis tends 
to be larger than if only one reference state is used. Thus, we have used a different set of importance-measure grid 
points, k = {7.00,6.00,5.00,4.00,3.75,3.50,3.25,3.00,2.75,2.50,2.25,2.00} x 10" 5 . This range is a bit smaller than 
the one we previously used for a single reference state and only extends to k = 2.0 x 10~ 5 instead of k = 1.0 x 10 -5 . 

In Fig. 11, we plot the relative difference to the NCSM result (horizontal line), for various extrapolating techniques, 
■Eg^o (left) or Eq 1 ^^ (right), as a function of increasing N max (hfl = 16 MeV). This is to be compared with Fig. 
4. The overall trend is the same between the two plots, indicating that the difference between IT-NCSM and NCSM 
calculations does not increase for the excited states, and that the difference is, in general, the same size as before 
(about 100 keV for N max = 14). In other words, the excitation spectrum can be calculated with the same degree 
of accuracy as for the ground-state. The extrapolation to N max = oo has a larger uncertainty than what has been 
shown in the earlier figures. The uncertainty is especially large for the extrapolation that uses the first-order energies 
(Eq 1 ^^). Considering only the first-order extrapolation, one sees that the predictions of the cubic- or quartic- 
polynomial extrapolation cither under- or over-estimate the NCSM result as iV max increases. Thus, the result is that 
the ensemble-averaging procedure at N max = oo has a large uncertainty since the two polynomials essentially force 
the N max = oo extrapolation in opposite directions (compared to the NCSM). 

An interesting comparison to make between using one or several reference states is to determine the behavior of 
the ground-state energy as a function of the number of basis states kept. This analysis is shown in Fig. 12, for 
Ml = 12 MeV. We chose that particular value of hfl, because the largest number of basis states are kept for this case. 
From the figure, one can deduce two interesting points. The first is that the additional reference states are selecting 
basis states that were previously not selected, as can be seen at the start of each N max space. These additional 
basis states tend to make the functional dependence of the ground-state energy as a function of the size of the basis 
approximately constant. It can also be seen that when multiple reference states are used, fewer states are needed 
than before in order to achieve the same ground-state energy as with a single reference state. The second point 
has to do with the lowering of the ground-state energy as a function of N max . Note that higher N max contributions 
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Figure 11. (color online) The figure shows the relative difference to the NCSM result (horizontal line), for various extrapolating 
techniques, Eg K=0 (left) or Eqx=o (right), as a function of increasing JV mal . In this case three reference states were used 
(J 71 = 1 + , 3 + , + ). We calculate the ground-state and the first two excited states. Note that the overall trend is the same for all 
three states, indicating that IT-NCSM performs equally well for excited states as for the ground-state. We have also performed 
the extrapolation to iV max = oo (black lines) by using the ensemble-averaging technique described in Sec. V. (Ml — 16 MeV). 



significantly lower the ground-state energy, when compared to simply adding more states in a single N max space. In 
other words, note that the drop in energy between N ma _ x = 12 — > 14 is larger than the drop in energy in just the 
Njoax = 12 space, which results from adding all the basis states that are kept. Such a feature could hold promise for 
doing configuration-interaction calculations, in which one- and two particle-hole excitations are created on top of a 
Hartree-Fock state determined in a small iV max space. Most of the binding energy is gained from adding the most 
significant configurations found in larger iV max spaces. However, we did notice in our IT-NCSM calculations that as 
Amax increases, previously discarded basis states in the lower -/V max spaces do become relevant at some stage and are 
added back into the basis by our basis evaluation procedure. This also explains why we typically use all states up to 
and including N niax = 4; those basis states will be added to the basis in any case by IT-NCSM. In fact, by the time 
we have completed our 7V max = 14 calculation, almost all of the N max = 6 states have been added to the basis, even 
though initially a fair number of those states were discarded at the start of the calculation. Such behavior makes 
sense, since we expect low-lying states to have components mostly found in the lower oscillator shells. 

We expect many of the same results to hold for multiple reference states, as shown for the single reference state 
calculations. In particular, we have determined that the Ml dependence, as shown in Fig. 8, still persists, when 
multiple reference states are used. This once again leads us to the conclusion that the dependence stems from k being 
inversely proportional to HQ. 



VII. CONCLUSION 



We have presented a detailed investigation of IT-NCSM calculations for 6 Li in which we have studied the dependence 
of IT-NCSM on various parameters. These include the behavior of IT-NCSM as a function of the model space N max , 
the HO energy M7, the extrapolating functions used for the two types of data sets (E^l or Eq£ ), the SRC 
momentum-decoupling scale A, as well as the influence on the basis selection procedure, when multiple reference 
states are used. The IT-NCSM calculations were then compared to NCSM calculations, as a way to estimate the 
efficiency of the procedure. We find that the extrapolations used in IT-NCSM, using cither the first- or second-order 
results (Eq^ =0 or ^Ore=o)i S^ ve similar results, even when different extrapolating functions are used. At N max = 14 
we find that the IT-NCSM extrapolated ground-state energies differ from the NCSM ground-state energies by about 
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Figure 12. (color online) The figure shows the ground-state energy at HQ — 12 MeV, as a function of the Logarithm (base 10) 
of the number of basis states kept. The curve with (+) signs corresponds to when only the ground-state is used as a reference 
state, wheres the curve marked (x) shows the behavior, when the lowest three states are used as reference states. Note that 
these energies correspond only to the first-order results, E^},, and that the two sets of kappa grid points are not identical. The 
figure shows that basis states coming from higher iV max spaces significantly lower the ground-state energy (note the decrease in 
the ground-state energy between each iV max space) as opposed to basis states in the same -/V max space that have a smaller value 
of re. Furthermore, using several reference states leads to an approximately constant dependence for the ground-state energy 
on the number of basis states kept. 



100-150 keV, whereas in the extrapolation to N max = oo the difference is about 250 keV. The IT-NCSM calculations 
show no quantitative difference for the two SRC momentum-decoupling scales that we used (A = 2.02 fm -1 and 
A = 1.50 fm -1 ) in terms of the extrapolation uncertainties or in regard to the difference to the full NCSM calculations. 
Two new features were seen in these calculations that have not been reported before: 1.) IT-NCSM calculations seem 
to deteriorate in quality as Ml increases, when the same re— grid is used (see Fig. 4 and Fig. 9); 2.) using several 
reference states leads to a better basis selection for the ground-state energy, than using just a single reference state 
(sec Fig. 12). We propose that future IT-NCSM calculations should use multiple reference states to select the basis 
states, use smaller re threshold limits as Ml increases and provide a reasonable uncertainty estimate of the iV max = oo 
extrapolated energies. 

In this work, we have tried to provide an improved analysis of uncertainties in IT-NCSM calculations. However, we 
must point out that the only way we are able to be confident in our IT-NCSM calculations is by direct comparison 
to full NCSM calculations. In situations where the NCSM calculations can be performed, IT-NCSM energy-spectra 
have compared favorably to NCSM [30, 31, 55]. However, in cases where no NCSM calculations exist, IT-NCSM 
calculations should be treated carefully, especially when ground-state energies are extrapolated to N max = oo. Recall 
that in order to extrapolate the ground-state energies to N max = oo two extrapolations must be performed. The 
first extrapolation is performed in a given N max space to estimate the NCSM result (re = 0). As N max increases 
the extrapolated ground-state energy typically underestimates the full NCSM ground-state energy, especially as Ml 
increases for a fixed re— grid. The underestimation of the ground-state energy in larger N max spaces inadvertently drives 
the extrapolation to -/V max = oo to a smaller binding energy than what is predicted with the full NCSM. Furthermore, 
one is free to choose a low-order polynomial function (or asymptotically-correct function) to extrapolate the ground- 
state energies in each iV max space. Each functional choice leads to a different prediction and different uncertainty for 
the ground-state energy. In our N max = oo extrapolation, we have taken into account the uncertainties of each N max 
ground-state energy and have determined the uncertainty at N max — oo over a range of Ml to be on the order of 
100-200 keV for A = 2.02 fm" 1 and on the order of 50-100 keV for A = 1.5 fm" 1 (see Fig. 8 and Fig. 9, respectively). 
The smaller uncertainty for A = 1.5 fm -1 is due to the softer nature of the underlying NN interaction. 
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